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Abstract. We study the statistical mechanics of a finite-dimensional non-linear Hamiltonian system (a chain of 
anharmonic oscillators) coupled to two heat baths (described by wave equations). Assuming that the initial conditions 
of the heat baths are distributed according to the Gibbs measures at two different temperatures we study the dynamics of 
the oscillators. Under suitable assumptions on the potential and on the coupling between the chain and the heat baths, 
we prove the existence of an invariant measure for any temperature difference, i.e., we prove the existence of steady 
states. Furthermore, if the temperature difference is sufficiently small, we prove that the invariant measure is unique and 
mixing. In particular, we develop new techniques for proving the existence of invariant measures for random processes 
on a non-compact phase space. These techniques are based on an extension of the commutator method of Hormander 
used in the study of hypoelliptic differential operators. 
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1. Introduction 

In this paper, we consider the non-equilibrium statistical mechanics of a finite-dimensional non- 
linear Hamiltonian system coupled to two infinite heat baths which are at different temperatures. 
We show that under certain conditions on the initial data the system goes to a unique non- 
equilibrium steady state. 

To put this new result into perspective, we situate it among other results in equilibrium 
and non-equilibrium statistical mechanics. First of all, for the case of only one heat bath one 
expects of course "return to equilibrium." This problem has a long history, and a proof of return 
to equilibrium under quite general conditions on the non-linear small system and its coupling to 
the heat bath has been recently obtained in [JP1-4] .Viewed from context of our present problem, 
the main simplifying feature of the one-bath problem is that the final state can be guessed, a 
priori, to be the familiar Boltzmann distribution. 

For the case of two heat baths, there are no results of such generality available, among other 
things precisely because one cannot guess in general what the steady state is going to be. Since 
we are dealing with systems on a non-compact phase space and without energy conservation, 
there is nothing like an SRB Ansatz for our problem [GC]. Worse, even the existence of any 
stationary state is not obvious at all. The only notable exceptions are problems where the small 
system and its coupling to the heat baths are linear. Then the problem can be formulated in terms 
of Gaussian measures, and approach to a steady state has been proved in this case in [RLL], 
[CL] and [OL] for Markovian heat baths and in [SL] for the general case. 

Our approach in the present paper will consist in using the spirit of [FKM] and [FK] to give 
a microscopic derivation of the equations of motion: under suitable assumptions, we will reduce 
the study of the dynamics of the coupled system (an infinite dimensional Hamiltonian system) 
to the study of a random finite dimensional dynamical system. However, we will not achieve 
the generality of [JP1-4]. Each heat bath is an infinite dimensional linear Hamiltonian system, 
in our case it will be chosen as the classical field theory associated with the wave equation. The 
small system is a non-linear Hamiltonian system with an arbitrary (but finite) number of degrees 
of freedom, in our case it is chosen as a chain of anharmonic oscillators with nearest neighbor 
couplings. The potential must be of quadratic type near infinite energies. The two heat baths 
are coupled respectively to the first and the last particle of the chain. The initial conditions of 
the heat baths will be distributed according to thermal equilibrium at inverse temperatures /3 L , 
(3 R . Integrating the variables of the heat baths leads to a system of random integro-differential 
equations: the generalized Langevin equations. They differ from the Newton equations of 
motion by the addition of two kinds of force, on one hand there is a (random) force exerted by 
the heat baths on the chain of oscillators and on the other hand there is a dissipative force with 
memory which describe the genuine retro-action from the heat bath on the small system. We will 
choose the couplings between the baths and the chain such that the random forces exerted by the 
baths have an exponentially decaying covariance. With this assumption (see [Tr]), the resulting 
equations are quasi-Markovian. By this, we mean that one can introduce a finite number of 
auxiliary variables in such a way that the evolution of the chain, together with these variables, 
is described by a system of Markovian stochastic differential equations. 

With this set-up, we are led to a classical problem in probability theory: the study of 
invariant measures for diffusion processes. For our problem, the main difficulties are as follows: 
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the phase space is not compact and the resulting diffusion process is degenerate and not self- 
adjoint. 1 The standard techniques used to prove the existence of invariant measures do not seem 
to work in our case and, in this paper, we develop new methods to solve this problem, which rely 
on methods of spectral analysis. Our proof of existence is based on a compactness argument, 
as often in the proof of existence of invariant measures. More precisely we will prove that the 
generator of the diffusion process, a second order differential operator, given in our problem, 
has compact resolvent, in a suitably chosen Hilbert space. This is done by generalizing the 
commutator method of Hormander, [Ho], used in the study of hypoelliptic operators. Similar 
methods have been used to study the spectrum of Schrodinger operators with magnetic fields, 
see [HM], [He], and [HN]. 

The restriction to a chain is mostly for convenience. Other geometries can be accommo- 
dated with our methods, and the number of heat baths is not restricted to two. Furthermore, the 
techniques developed in this paper can be applied to other interesting models of non-equilibrium 
statistical mechanics, for example, an electric field acting on a system of particles [R-B, in prepa- 
ration]. 

2. Description of the Model and Derivation of the Effective Equa- 
tions 

In this section we define a model of two heat baths coupled to a small system, and derive the 
stochastic equations which describe the time evolution of the small system. The heat baths 
are classical field theories associated with the wave equation, the small system is a chain of 
oscillators and the coupling between them is linear in the field. 

We begin the description of the model by defining the "small" system. It is a chain of 
d-dimensional anharmonic oscillators. The phase space of the chain is R 2dn with n and d 
arbitrary and its dynamics is described by a C°° Hamiltonian function of the form 

H s (p,q) = £^+£W + £^ 2) (<M + i) = E| + y (?), (2-1) 

j=\ j=l i=\ j=l 

whereg = (q t , . . . , q n ), p = (p u . . . ,p n ), with^, ^ G R d . 

The potential energy will be assumed "quadratic + bounded" in the following sense. We 
let T denote the space of C°° functions on R dn such that for all multi-indices a and all F G T , 
the quantity d a F(q) is bounded uniformly in q G H dn . Then our hypotheses are 

HI) Behavior at infinity: We assume that V is of the form 

V(q) = ±(q-a,Q(q-a))+F(q), 

where Q is a positive definite (dn x dn) matrix, a is a vector, and d ( V )F G T for 
i = 1 , . . . , n and v = 1 , . . . , d. 

1 The diffusion is non-selfadjoint because the diffusion process is not time-reversal invariant. But the generator L 
satisfies L*T = TL where T changes the sign of the momenta. 
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H2) Coupling: Each of the (d x d) matrices 

M iji+1 (q) = V ft V ft+1 C^ 2) ( ?i , ?4+1 ), i = l,...,n-l, 

is either uniformly positive or negative definite. 

Remark. The first hypothesis makes sure the particles do not "fly away." The second hypothesis 
makes sure that the nearest neighbor interaction can transmit energy. As such, this condition is 
of the hypoelliptic type. 

Example. A typical case (in dimension d) covered by these hypotheses is given by 

Uj(q) = Q 2 + 5sin( v / l + ? 2 ) , tfW) = {q - q'f + ^\ + {q - q') 2 ) . 

As a model of a heat bath we consider the classical field theory associated with the d- 
dimensional wave equation. The field cp and its conjugate momentum field tt are elements of 
the real Hilbert space H = (R d ) © (R d ) which is the completion of C °° (R d ) © C °° (R d ) 
with respect to the norm defined by the scalar product: 

((i)>(t)) H = /dx(|V^)| 2 +Kx)| 2 ) . (2.2) 
The Hamilton function of the free heat bath is 



H^) = \ j dx (|V^(x)| 2 +|7r(x)| 2 ) , 



and the corresponding equation of motion is the ordinary wave equation which we write in the 
form 



where 

L = [A 

Let us turn to the coupling between the chain and the heat baths. The baths will be called 
"L" and "R", the left bath couples to the coordinate q 1 and the right bath couples to the other 
end of the chain (q n ). Since we consider two heat baths, the phase space of the coupled system, 
for finite energy configurations, is R 2dn x H x H and its Hamiltonian will be chosen as 

H(p, q, <p L , tt l , v? R , tt r ) = H s (p, q) + H B (<p L , tt l ) + H B (<p R , tt r ) 

f f ( 2 -3) 

+ Qi- J d *Ph( x )V<P L (x) +q n ■ J dxp R (x)Vip R (x) . 

Here, the Pj(x) G L 1 (R d ) are charge densities which we assume for simplicity to be spherically 
symmetric functions. The choice of the Hamiltonian Eq.(2.3) is motivated by the dipole 
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approximation of classical electrodynamics. For notational purposes we use in the sequel the 
shorthand 

We set = (c^,...,^)^ G {L,R}, with 

a<T\k)B 

Here and in the sequel the "hat" means the Fourier transform 



o 



—ik-x 



With this notation the Hamiltonian becomes 

H(p, q, </> L , R ) = H s (p, q) + H B ((j) L ) + H B (0 R ) + q x ■ (0 L , a L ) H + q n ■ (0 R , a R ) H , 
where H B (<fi) = \ ||0||^. We next study the equations of motions. They take the form 

Qj(t) = Pj(t) , j = l,...,n, 
Pi(t) = -V^(g(*))-(0L(*),a L ) w , 

Kfaftn, 7=2,...,n-l, 

(2.4) 



= - V «^(9(*)) ' j = 2,...,n-l, 



P„(*) = -V^(?(*))-(0r(*W)«, 
L (t) = £(0 L (t) + a L •?!(£)) , 

=£Wt) + a R -g B (t)) • 

The last two equations of (2.4) are easily integrated and lead to 

Mt) = e £ V L (0)+ fdsLe L ^a^-q x {s), 
Jo 

Mt) = ^Vr(0)+ fdsU^'^^-qJs), 
Jo 

where the ^(0), % E {L, R}, are the initial conditions of the heat baths. Inserting into the first 
2n equations of (2.4) gives the following system of integro-differential equations 



Qj(t) = Pj(t) , j = l,...,n, 

Pi(t) = -^ qi V(q(t))-(M0),e- Ct a L ) n - f* ds D L (t - s)q l (s) , 

Jo 

P^) = -V gi ^(g(*)) , j = 2,...,n-l, 

Pn(t) = -V q V(q(t))-(M0U-''W) n - [ ds D R (t - s)q n (s) 

Jo 



(2.5) 
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where the d x d dissipation matrices D\ (t — s), i £ {L, R}, are given by 
D^ v \t-s) = (a^ , -Ge^'-^) 

= ~\ 5 ^j dk \Mk)\ 2 \k\sm(\k\(t-s)). 
The last expression is obtained by observing that 



e 



u = ( cos(|fc|t) \k\~ l sin(\k\t)\ 
-\k\sm(\k\t) cos(]/c|£) / 



written in Fourier space. 

So far we only discussed the finite energy configurations of the heat baths. We will assume 
that the two reservoirs are in thermal equilibrium at inverse temperatures (3 L and (3 R . This means 
that the initial conditions $(0) = {0 L (O), R (O)} are distributed according to the Gaussian 
measure with mean zero and covariance (<j>i(f)<j>j(g)) = 8 i j(l/P i )(f,g) n . (Recall that the 
Hamiltonian of the heat baths is given by J2ie{h r}(^> ^ we assume that the coupling 

functions arein7i, i G {L, R}, and v e {1, • • • , d} then the ^(t) = <p i (0)(e~ iit a i ) become 
(i-dimensional Gaussian random processes with mean zero and covariance 

&(*)£;(«)> = koJ C S ~ s ) > G {L, R} , (2.6) 
and the d x d covariance matrices C^t — s) are given by 

C^\t-s) = (ofyNoj^ 

= I^ydk^fc^cosflfcKt-a)). 

The relation 

C 4 (f) = !>,(*), (2.7) 

which is checked easily by inspection, is known as the fluctuation dissipation theorem. It is 
characteristic of the Hamiltonian nature of the system. After these assumptions and transforma- 
tions, the equations of motion (2.5) become a system of random integro-differential equations 
on R 2dn which we will analyze in the sequel. 

Finally, we impose a condition on the random force exerted by the heat baths on the chain. 
We assume that 

H3) The covariances of the random processes £^(£) with i e {L, R} satisfy Cf'^\t — s) = 

<W Em=l t^'""" 1 ' With %m > ^d \ m > 0. 

This can be achieved by a suitable choice of the coupling functions p^x), for example 

M 



m) =<*»*. n ( F + i,j./ 2 - 
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where all the r y i m are distinct. To keep the notation from still further accumulating, we choose 
M the same on the left and the right. We will call the random process given by Eq.(2.5) 
quasi-Markovian if Condition H3 is satisfied. Indeed, using Condition H3 together with 
the fluctuation-dissipation relation (2.7) and enlarging the phase space one may eliminate the 
memory terms (both deterministic and random) of the equations of motion (2.5) and rewrite 
them as a system of Markovian stochastic differential equations. 

By Condition H3 we can rewrite the stochastic processes £-(t) as Ito stochastic integrals 



m=l V J-oo 



where the w i m (s) are (/-dimensional Wiener processes with covariance 



E 



" i^Ut') - rfUsj)] = \A,» 6 m,m>M n [S',t']\ , (2.8) 



where s < t and s < t' , E is the expectation on the probability space of the Wiener process and 
| • | denotes the Lebesgue measure. We introduce new "effective" variables r L m , r R m G R d , 
with m = 1 , . . . , M, which describe both the retro-action of the heat bath onto the system and 
the random force exerted by the heat baths: 

r L , m (t) = \l, m l^ m f^se-^- s \,{s) -A L>m yS^ f ^ e^^dv^s) , 

r R , m (t) = ^,mlR,mJ*dse-^-^ qn (s) - \ R , m ^^ f J-^ m{t - s) ^ m {s). 

We get the following system of Markovian stochastic differential equations 

= Pj(t)dt, j = l,...,n, 

M 
m=l 

d Pj (t) = -V q V(q(t))dt, j = 2,...,n-l, 

M 

dP n (t) = -V q V(q(t))dt + £ r Rtm (t)dt , (2.9) 

m=l 



d^L,m(*) = -7 L ,m^L,mW d ^ + ^L,m7 L ,m9lWdt- A^ m Ji^d«; L m (t) , 



dr R,mW = -7R, m r R,m(^) di + A R,m7R, m ?n(^ ~ V,ml/^ d ™R,m(*) > 
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which defines a Markov diffusion process on R ^ ' . This system of equations is our main 
object of study. Our main results are the following: 

Theorem 2.1. If Conditions H1-H2 hold, there is a constant A* > 0, such that for |A L m \, 
| A R m \ E (0, A*) with m = 1, . . . , M, the solution of Eq.(2.9) is a Markov process which has 
an absolutely continuous invariant measure fi with a C°° density. 

Remark. In Proposition 3.6 we will show even a little more. Let h (f3) be the Gibbs distribution 
for our system when the heat baths are both at temperature 1 / f3. If h denotes the density of the 
invariant measure found in Theorem 2.1, we find that h/h (/3) is in the Schwartz space S for 
all (3 < min(/3 L , f3 R ). This mathematical statement reflects the intuitively obvious fact that the 
chain can not get hotter than either of the baths. 

Concerning the uniqueness and the ergodic properties of the invariant measure, our results 
are restricted to small temperature differences. We have the following result. 

Theorem 2.2. If Conditions H1-H2 hold, there are constants A* > and e > such that 

for l\,ml> W,m\ e (°» A *) with 771 = 1, • • • , M, and \f3 L - P R \/(P h + P R ) < e, the Markov 
process (2.9) has a unique invariant measure and this measure is mixing. 

Remark. The restriction on the couplings between the small system and the baths A L m , A R m 
is non-perturbative: it is a condition of stability of the coupled system small system plus heat 
baths. Indeed, the baths have the effect of renormalizing the deterministic potential seen by the 
small system. The constant A* depends only on the potential V(q): if the coupling constants 
A L m , A R m are too large, the effective potential ceases to be stable and, at least at equilibrium 
(i.e., for P h = f3 R ), there is no invariant probability measure for the Markov process (2.9), 
but only a a-finite invariant measure (see Eq.(3.7) and Eq.(3.9)). This restriction is related to 
Condition HI on the potential: for potentials which grow at infinity faster than quadratically, 
this restriction would not be present (see [JP1]). On the other hand, the restriction on the 
temperature difference is of perturbative origin. 

A more physical formulation of the results of Theorem 2.2 is obtained by going back to 
the Eq.(2.5), which expresses all the quantities in terms of the phase space of the small system 
and the initial conditions $(0) of the heat baths. Let us introduce some notation: For given 
initial conditions $(0), we let Q t $( )(p, q) denote the solution of Eq.(2.5). Finally, define 

i/(dp, dq) = / //(dp, dq, dr) , 

JrE'R 2dM 

where n is the invariant measure of Theorem 2.1. 

Corollary 2.3. Under the hypotheses of Theorem 2.2, the system Eq.(2.5) reaches a stationary 
state and is mixing in the following sense: For any observables F,G E L 2 (R 2dn , i/(dp, dq)) and 
for any probability measure ^ (clp, dq) which is absolutely continuous with respect to z/(dp, dq) 
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we have 




lim / I /(dp,dq)((Foe t)$o )(p,g)G(p,g)> = /i/(dp, dq)F(p, q) h(dp, dq)G(p, g) . 



Here, (•) denotes the integration over the Gaussian measures of the two heat baths, introduced 



We explain next the strategy to prove these results. Our proof is based on a detailed study 
of Eq.(2.9). Let x = (p, q, r). For a Markov process x(t) with phase space X and an invariant 
measure jii(dx), its ergodic properties may be deduced from the study of the associated semi- 
group T* on the Hilbert space L (X, /i(dx)). To prove the existence of the invariant measure in 
Theorem 2. 1 we proceed as follows: We consider first the semi-group T* on the auxiliary Hilbert 
space H Q = L 2 (X, / u (dx)), where the reference measure £i (dx) is a generalized Gibbs state 
for a suitably chosen reference temperature. Our main technical result consists in proving that 
the generator L of the semi-group T* on 7i and its adjoint L* have compact resolvent. This is 
proved by generalizing the commutator method developed by Hormander to study hypoelliptic 
operators. From this follows the existence of a solution to the eigenvalue equation (T 1 )* g = g 
in H and this implies immediately the existence of an invariant measure. To prove Theorem 2.2 
we use a perturbation argument, indeed at equilibrium (i. e. , for f3 L = f3 R ) the invariant measure is 
unique and is a simple eigenvalue of the generator L in H . Using the compactness properties 
of L, we show that is a simple eigenvalue of the generator L in H for \(3 L — (3 R \/((3 L + /? R ) 
small enough. And this can be used to prove the uniqueness claim of Theorem 2.2, while the 
mixing properties will be shown by extending the method of [Tr] . 

This paper is organized as follows: In Section 3 we prove Theorem 2.1 and Theorem 2.2 
except for our main estimates Proposition 3.4 and Proposition 3.5 which are proven in Section 
4. In Appendices A, B and C, we prove some auxiliary results. 

3. Invariant Measure: Existence and Ergodic Properties 

In this section, our main aim is to prove Theorem 2.1 and Theorem 2.2. We first prove 
some basic consequences of our Assumptions HI and H2. In particular, we define the semi- 
group T* describing the solutions of Eq.(2.9) on the auxiliary Hilbert space H described in 
the introduction. Furthermore we recall some basic facts on hypoelliptic differential operators. 
Once these preliminaries are in place, we can attack the proof of Theorem 2. 1 and Theorem 2.2 
proper. 



(2.10) 



earlier. 
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3.1. Existence and fundamental properties of the dynamics 

Let X = rM™+ m ) an d wr ite the stochastic differential equation (2.9) in the abbreviated form 

dx(t) = b(x(t)) + adw(t) , (3.1) 

where 

(i) b is a C°° vector field which satisfies, by Condition HI, 

sup \d a b(x)\ < oo , 

xex 

for any multi-index a such that \a\ > 1. In particular 

fl= Hdiv&L, < oo. (3.2) 

(ii) a : K 2dM - X is a linear map. We also define 

D = X -ao T > 0. (3.3) 

(iii) w G W = C(R; R 2dM ) is a standard 2<iM-dimensional Wiener process. 
The Eq.(3.1) is rewritten more precisely as the integral equation 

£(t,w;x) = x + [ dsb(£(s,w,x)) + a(w(t)-w(0)) . (3.4) 
Jo 

It follows from an elementary contraction argument (see e.g. [Ne], Thm 8.1) that (3.4) has a 
unique solution 

R 9 t h-> x (t) = £(t, w; x) e C(R; X) , 

for arbitrary initial condition x E X and w G W. 

The difference u;(£) — ty(0) has the statistics of a standard Brownian motion and we 
denote by E[-] the corresponding expectation. By well-known results on stochastic differential 
equations, this induces on £ (t , w; x) the statistics of a Markovian diffusion process with generator 

V ■ DV + b(x) - V . (3.5) 

More precisely (see [Ne] Theorem 8.1): Let C OQ (X) denote the continuous functions which 
vanish at infinity with the sup-norm and let T l be the a-field generated by x and {w(s) — 
w(0) ;0<s<t}, then for < s < t and / G C^X) we have 

E[/(x(t))|^ s ] = T^fixis)) a.s., (3.6) 

where T* is a strongly continuous contraction semi-group of positivity preserving operators on 
C 00 (X) whose generator reduces to (3.5) on C^°(X). 
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In the sequel we denote by L the differential operator (3.5) with domain D(L) = C£°(X). 

To prove the existence of an invariant measure we will study the semi-group T* or rather 
an extension of it on the auxiliary weighted Hilbert space H described in the introduction. To 
define H Q precisely, we consider the "effective Hamiltonian" 



M ( 1 r tm 1 r Rm \ 

G(p, q, r) = H s (p, q) + J2\ T2 + 72 ~ 9i ' r L,m ~ Q n ' r R, m • ( 3 - 7 ) 

m=1 \ A L,m A A R,m L J 

We note that, due to Condition HI, G(x) — > +oo as \x\ — > oo as long as |A L m \, |A R m | < A* 
for some A* depending only on the potential V(q). 

We choose further a "reference temperature" f3 Q , which is arbitrary subject to the condition 

p < 2min(/3 L ,/3 R ). (3.8) 

For example we could take (3 Q as the inverse of the mean temperature of the heat baths: 
[3q 1 = (/?£" 1 + (3 R ~ l )/2. For the time being, it will be convenient not to fix (3 . Then, we let 

H = L 2 (X,Z - 1 e _ftG dx) , (3.9) 

and we denote (-,-)n anc * II ' llw me corresponding scalar product and norm. 

Remark. With a proper choice of Z , it is easy to check that the quantity ZQ 1 e~ /3 ° G ^ q,p ' r ' dx is 
the invariant measure for the Markov process Eq.(2.9) when f3 L = f3 R = /3 Q and | A L m | , | A R m \ < 
A*. " ' 

Lemma 3.1. If the potential V satisfies Condition HI and if (3 < 2min( / 3 L , (3 R ) there is a 
A* > such that if the couplings satisfy |A Lm |, |A R m \ e (0, A*), then the semi-group T* given 
by Eq.(3.6) extends to a strongly continuous quasi bounded semi-group T^ on H : 



\ r r t II < p at 

■J- 1-/ c 

I /TO II /TO — 



where a is given by 



a = dJ2[ 7 L ,m 

m=l V V 



1 ^/(P L -P /2)P /2 \ ( 1 y/(0 R -0 o /2)0 o /2 \ 

2 — k — ) + ^{2 — s — j ; - ° 

(3.10) 

The generator L nQ ofT^ Q is the closure of the differential operator L with domain given by 

M ,2 

m=\ L 
M ,2 

+ E Jb f^ - /3rWr -) • v — (3 - n) 

m=l ^ R 
M M 

+ E r L,m • V pi + L s + r R,m ' ^ Pn , 

m=l m=l 
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with the abbreviations 

W Um = ^L,m r L,m ~ ?1 > W R,m = A R 2 m r R,m ~ Qn , ( 3 - 12 ) 

and where L s is the Liouville operator associated with the Hamiltonian H s (q, p): 

n 

Moreover, T^ o is positivity preserving: 

T'nJ > if / > 0, (3.14) 

and 

2^1 = 1. (3-15) 
Remark. We have a = if only if /3 L = /3 R = /? . 

Proof. The proof uses standard tools of stochastic analysis and is given in Appendix A. 

Having shown a priori bounds using Condition HI, we will state one basic consequence 
of Condition H2. We recall that a differential operator P is called hypoelliptic if 

sing supp u = sing supp Pu for all u E V'(X) . 

Here V'(X) is the usual space of distributions on the infinitely differentiable functions with 
compact support and for u € V'(X), sing supp u is the set of points x £ X such that there is 
no open neighborhood of x to which the restriction of u is a C°° function. 
Let P be of the form 

J 

P = E^ + y o, (3-16) 

3 = 1 

where Y-,j E {0, . . . , J} are real C°° vector fields. Then by Hormander's Theorem, [Ho], Thm 
22.2.1, if the Lie algebra generated by Y-, j e {0, . . . , J} has rank dim X at every point, then 
P is hypoelliptic. 

Differential operators arising from diffusion problems are of the form (3.16). Let L be the 

rp 

differential operators given in Eq.(3.1 1), let L denote its formal adjoint, then one may easily 
check that Condition H2 implies that any of the following operators 

L, L T , d t + L, d t + L T , 

satisfies the condition of Hormander's Theorem and thus is hypoelliptic. As an immediate 
consequence we have: 
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Corollary 3.2. If Condition H2 is satisfied then the eigenvectors of L and L T are C°° 
functions. 

Next, letP(t, x, E), t > 0, x G X,E G B denote the transition probabilities of the Markov 
process £(£, w; x) solving the stochastic differential equation (2.9) with initial condition x, i.e., 

P(t,x,E) = P(£{t,w:x)eE) , 

where P denotes the probability associated with the Wiener process. Then by the forward and 
backward Kolmogorov equations we obtain 

Corollary 3.3. If Conditions HI and H2 are satisfied then the transition probabilities of the 
Markov Process £(t,w; x) have a smooth density 

P(t, x, y) G C°°((0, oo)xX xX) . 



3.2. Proof of Theorem 2.1 and Theorem 2.2 

After these preliminaries we now turn to the study of spectral properties of the generator L Ho 
of the semi-group TL . 

The proof of the existence of the invariant measure will be a consequence of the following 
key property which we prove in Section 4. 

Proposition 3.4. If the potential V satisfies Conditions HI, H2 and if (3 < 2min( / 3 L , /3 R ) 
there is a A* > such that if the couplings satisfy |A L m |, |A Rm | G (0, A*) then both L Ho and 
L^ o have compact resolvent. 

A useful by-product of the proof Proposition 3.4 are some additional smoothness and 
decay properties of the eigenvalues of L Hq and L^ o on H . 

Proposition 3.5. Let g denote an eigenvector of L Hq or L^. If the assumptions of 
Proposition 3.4 are satisfied then we have 

ge -PoG/2 e ^ 

where S(X) denotes the Schwartz space. 

Using these results, we come back to the Markov process defined by the Eqs.(2.9), and 
whose semi-group T l was defined in Eq.(3.6). We prove the existence of an invariant measure 
with a smooth density and give a bound which shows that, in some sense, the chain does not get 
hotter than the hottest heat bath. 

Proposition 3.6. Under the assumptions H1-H2 there is a X* > such that if the couplings 
satisfy |A L \, |A R | G (0, A*) the Markov process T l has an invariant measure /j, which is 
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absolutely continuous with respect to the Lebesgue measure. Its density h satisfies the following: 

h&xp(/3G) e S(X) for all (3 < min(/3 L , /3 R ). 

Proof. The function 1 is obviously a solution of Lf = with L defined in Eq.(3.11). Note 
next that the function 1 is in H , as is seen from Eq.(3.9) (if |A L m \ and |A R m \ are sufficiently 
small). Since, by Proposition 3.4, the operator L n has compact resolvent on Ti , it follows 
that is also an eigenvalue of L* Hq . Let us denote the corresponding eigenvector by g. We will 
choose the normalization (g, l) w = 1. We assume first that g > 0. Then the function 

h(x) = Z^g(x)e-^ , (3.17) 

with P Q and G defined in Eqs.(3.8) and (3.7), is the density of an invariant measure for the process 
T l : Indeed, we note first that |H| L i(x,dx) = (l)#)w * s ^ n ^ e an ^ thus /i(dx) is a probability 
measure. 

Let E be some Borel set. Then the characteristic function \e °f & belongs to H . We 

have 

J ^dx)T\ E = V/dx e-* G Mg(z)T t XE 
= KE) , 

and therefore //(dx) is an invariant measure for the Markov process (2.9). 

To complete the first part of the proof of Proposition 3 .6 it remains to show that g > 0. We 
will do this by checking that h > 0. We need some notation. Let L T denote the formal adjoint 
of L. Then one has L J h = 0. This follows from the identities 

JdxfL T h = V JdxfL T {ge-^ G ) 

= V J dx {Lf)ge-^ G 

= {Lf,9) Ho = (f^n g) HQ = (f,0) Ho =0, 

which hold for all / e C£°(X). Consider now the semi-group T* acting on the space C 0Q (X) 
defined at the beginning of Section 3. Since T is a Markovian semi-group, we haveT*l = 1 
and T* is positivity preserving. The operator T* induces an action (T*)* defined on the dual 
space (X) which consists of finite measures. Since T* is Markovian, (T*)* maps probability 
measures to probability measures. Furthermore, if a measure v has a density / in L*(X, dx), 
then (T*)*i/ is a measure which has again a density in L (X, dx): Indeed, by Corollary 3.3 the 
transition probabilities of the Markov process P(t, x, y) are in C°° ((0, oo) x X x X) . If we 
denote by (T*) T the induced action of (T*)* on the densities, we have for g > 0, 



(T t y(g(x)dx) = J dyg(y)P(t,y,dx) = dx J dy g(y)P(t, y, x) = ((T'fg) (x)dx , 
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and IK^Hl^IMIl.. 

Coming back to the invariant density h, we know that 

(T'fh = h . 

We next show (T t ) T |/^| = \h\. Since \h\±h> 0, we have (T*) T (|/i| ± h) > 0. This can be 
rewritten as 

|(T*) T /i| < (T'flh] . 

Therefore, 

\h\ = \{T l fh\ < (T*) T |/i| . 

Since {T 1 ) 1 preserves the L -norm, we conclude that 

\h\ = (T^l/il . (3.18) 

This shows the existence of an invariant measure. 

Now, by Proposition 3.5, we have h exp(/3G/2) G S(X) for all (3 < 2 min(,9 L , (3 R ) and so 
for/3 < min(/3 L , f3 R ) it follows that hexp(PG) G «S(X). This concludes the proof of Proposition 
3.6. 

We next prove the uniqueness of the invariant measure and the ergodic properties of the 
Markov process. We start by fixing an inverse temperature (3 Q . If (3 L = (3 R = (3 , the two 
heat baths are at the same temperature, and the equilibrium state of the system is known, since 
it is given by the generalized Gibbs distribution ZQ 1 e~ ,3oG . For the equilibrium case, this 
distribution is the unique invariant measure. The existence is obvious from what we showed 
for the case of arbitrary temperatures. To show uniqueness, assume that there is a second 

rp 

invariant measure. Since L is hypoelliptic, then by Corollary 3.2 this measure has a smooth 
density. Since different smooth invariant measures have mutually disjoint supports and e~^ G 
has support everywhere, uniqueness follows. If the invariant measure is unique, it is ergodic 
and hence, (see [Yo] and [Ho]) is a simple eigenvalue of L Hq . 

The case of different temperatures will be handled by a perturbation argument around the 
equilibrium situation we just described. This perturbation argument will take place in the fixed 
Hilbert space H defined in Eq.(3.9). Thus, we will consider values of f3 h and f3 R such that 

for some small e > (which does not depend on (3 ). 

We first show that remains a simple eigenvalue of the generator L Ha when the temperature 
difference satisfies (3.19) for a sufficiently small e. 

Lemma 3.7. Under the assumptions H1-H2 there are constants A* > and e > such that 
if the couplings satisfy |A L m |, |A R m \ G (0, A*) and moreover f3 h , f3 R satisfy (3.19), then is a 
simple eigenvalue of the generator L Ho . 

Proof. It will be convenient to work in the flat Hilbert space L 2 (X, dx). Note that fC = 

exp (— P G/2)Lexp (+/3 G/2) is a function /C = /C(/3 L , /3 R , f3 ). We write this as 

K(p L ,p R ,P ) = K(p ,p ,p ) + 8Z , 
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where 



One finds 



5 = 



M ,2 



m=l 



ft, 



M x2 
_j_ A R,m7R, m 



m=l 

M 



— V 

A) 



L,m 



R,m 



m=l 



m=l 



and 



M \2 



m=l 



M x 2 

A R,m7R, 



m=l 



,m * r R , m 1 * r Rjm ' ' R,m ^ 2 R > m 



Furthermore, by Proposition 3.4, i? = (1 — /C(/? , /9 , /3 )) _1 is a compact operator, and 
therefore the simple eigenvalue 1 of i? is isolated. From now on we assume for convenience 
that a = a(P L , (3 R ) is strictly smaller than one. Note that this is no restriction of generality: if 
a e [n — 1, n) with n > 1, we replace (1 — fC)~ l by (1 — ^/C) _1 in the following discussion. 

We show next that the resolvent R((3 L , [3 R , /3 ) = (1 — /C(/3 L , /3 R , @o))~ depends analyti- 
cally on the parameter 5. It is convenient to write the perturbation Z as 

JV 
3 = 1 



3 " 3 

and AT" = MM. With the matrix notation 



where the E- and F ? - are of the form const.dj,,) or const. W^, i G {L, R}, m = 1, . . . , M, 

x. notatic 



E T = (E l ,...,E N ) 



\FnJ 



we can write Z as Z = E F. We will use the following resolvent formula 

R(p L , p R , p ) = R (l + 5R E T (1 - R E T ^j ~ l FR^j . (3.20) 



Anharmonic Chains Coupled to Two Heat Baths 



17 



To justify Eq.(3.20) we have to show that for 6 small enough the operator- valued matrix 
(l — 5FR E T ) is invertible. It is enough to show that F-R Q E k is a bounded operator, for 
all j, k. For this we decompose 1 — /C(/? , /9 , f3 ) into its symmetric and antisymmetric parts: 

l-K(P ,P ,P ) = X + iY , 

where 

M i2 / O R 1 

x = , + E ^. 2^ + _„ 1 

m=l Z \ PO ^ L.» 

m=l 

From the simple estimates 



A) rR ' m 2 A| >m 



F,-/|| 2 < (/,*/), 



3 



which hold for / e C^°(X), and since X is a strictly positive operator we see that ly ^ 2 
and F^X -1 / 2 are bounded operators. From the identity 

E 3 R Q F k = Ej(X + iY)~ x F k 

= E 3 x- y i 1 (\ + ix- x i 1 Yx- x i 1 y x x- x i 1 F k , 

we see that E-R Q F k are bounded operators for all j, k. Therefore the r.h.s of Eq.(3.20) is well 
defined for sufficiently small 5. An immediate consequence of the resolvent formula (3.20) is 
that for sufficiently small 8 the spectrum of R(f3 L , f3 R7 f3 ) has the same form as the spectrum 
R : 1 is an eigenvalue and there is a spectral gap and, in particular 1 is a simple eigenvalue. 
This concludes the proof of Lemma 3.7. 

Next we use this lemma to prove uniqueness of the invariant measure. We have the 
following 

Theorem 3.8. Under the assumptions H1-H2 there are constants A* > ande > such that 
if the couplings satisfy |A L |, |A R | G (0, A*) and the temperatures satisfy |/5 R — /3 L |/(/3 L + 
(3 R ) < e, the Markov process T l has a unique (and hence ergodic) invariant measure. 

Proof. The proof uses a dynamical argument. By Proposition 3.4 we have in the Hilbert space 
Hq (with P given as in (3.19)) the eigenvalue equation T^ o l = 1 and (T^)*g = g. Let the 
eigenvectors be normalized such that (g,l) Ho = 1. By Lemma 3.7, is a simple eigenvalue of 
the generator L nQ if (/3 R — (3 h ) / (f3 L + f3 R ) is small enough and by Proposition 3.6, the measure 
jit(dx) = ZQ l gexp (—/3 G) is an invariant measure for the Markov process. It is absolutely 
continuous with respect to the Lebesgue measure which we denote by A. 

Assume now that v is another invariant probability measure. By the hypoellipticity of L 
it must have a smooth density. Therefore there is a Borel set A C X, which we may assume 
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bounded, with the following properties: we have u(A) > and A (A) > but (j,(A) = 0, 
because the measures have disjoint supports. Let Xa denote the characteristic function of the set 
A. By the pointwise ergodic theorem, see [Yo] and [Ho], we have, denoting a s (x) the ergodic 
average a s (x) = (1/s) f* dt T*x A (x), 

lim a Ax) = u(A) , i/-a.e. . (3.21) 

s^oo 

Since T l is a contraction semi-group on B{X,B), and < 1 we find HcrJI^ < 1, for all 
s > 0. From the easy bound 

\Ws\\ho ^ IkJIoo ' 

we see that the set {a s , s > 0} is a bounded subset of H and hence weakly sequentially 
precompact. Therefore there is a sequence s n | oo such that 

w-lim a s =cr* 

n — >oo ™ 

where w-lim denotes the weak limit in H . Since T u is a bounded operator for all u > 0, we 
have 

w-lim TV. =T u a* . 

We next show 

T u a* = a* . (3.22) 

Indeed, 

f s n +u 



T u a Sn (x) = - / diT^s) 

b n Ju 

- / dtT t XA (x) + — / d*T* XA (x). 

'n JO S n J s„ 



n Ju (3 23) 



S 



The last two terms in (3.23) are bounded by u/s n and we obtain T u a* = a* for all u > by 
taking the limit n — > oo in (3.23). 

Therefore, a* is in the eigenspace of the eigenvalue 1 of T-^ , t > and so, by Lemma 3.7 
a* = cl. To compute c we note that c = cr*)^ and, using the invariance of the measure, 
we get 

c = lim L_L f n dtT^ XA ) 

n ^°° \ S uJ0 J Ho 



l im L f n dt f fi (dx)T t x A = f i(A) 



So we have c = n(A) and /z(-A) = by hypothesis. Using this information, we consider 
(x A i a sJn - We have ' on one hand ' 

lim (xa^sJh = (Xa.Owo = °' 
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and on the other hand we have, by Eq.(3.21) and by the dominated convergence theorem, 

Jjim (Xa^sJho = J dxv(A)x A Z^e-P° G >0, 

and this is a contradiction. This shows that there is a unique invariant measure for the Markov 
process T* and as a consequence the measure is ergodic. 

We will now strengthen the last statement by showing that the invariant measure is in fact 
mixing. This will be done by extending the proof of return to equilibrium given in [Tr]. 

Proposition 3.9. Assume that the conditions of Theorem 3.8 are satisfied. Then the invariant 
measure /i(dx) for the Markov process T* is mixing, i.e., 

Urn J Mdx)/(*)TV(*) = J /x(dx)/(x) J ^dx)g(x) , 

forallf,geh 2 (X^(dx)). 

Proof. We denote H = L 2 (X, /i(dx)) and its scalar product by (-,•)■« an ^ by || • \\ n its norm. 
By [Yo], Chap XIII. 1, Thm 1, T* defines a contraction semi-group on H. Since T* is a strongly 
continuous semi-group on C OQ (X) (see [Ne]) and since C^X) is dense in H, we can extend 
T* to a strongly continuous semi-group on 7Y. The property of mixing is equivalent to 

w-limT^/ = (1, f) n for aU/eW. (3.24) 

By a simple density argument it is enough to show (3.24) for a dense subset of H. Let C 2 (X) 
denote the bounded continuous function whose first and second partial derivatives are bounded 
and continuous. Then— [GS], Part II, §9— if / G C 2 (X), then T l n f e C 2 (X) and for any 
t < oo, T^f is uniformly differentiable w.r.t. t G [0, r] and 

where L is the differential operator given in (3.11). Let / G C (X). Using the fact, see 
Proposition 3.5 and Proposition 3.6, that the density of the invariant measure is of the form 
h(x) = ge~P° G/2 with g G S(X), we may differentiate under the integral and integrate by parts 
and using the invariance of the measure we obtain 

±\Kf\\ 2 H = (LTy,Ty) n + ( T y,LTy) H = -J2 2 -^f^\\d r ^Ty\\ 2 H , 

;e{L,R} Pl 

m€{l,...,M} 

ve{\,...,d} 

(3.25) 
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where d i» is the differential operator with domain C 2 (X). Thus ||T^/||^ is decreasing, 
bounded below and continuous and so lim^^ ||T^/||^ exists. As a consequence we find 

\\d rM T^ff H E L\[0,oo),dt). (3.26) 

Following [Tr] and [Br], we call a sequence {t n } a (*)-sequence if t n | oo and 

lim \\d r(v) T%f\\ 2 n = Q. (3.27) 

71 »00 i,m 

The existence of (*)-sequences for our problem follows easily from (3.26). Further we define 
an almost (*)-sequence as a sequence s re f cxo for which there exists a (*)-sequence {t n } with 
\s n — t n \ — > as n — > oo. As in [Tr] we next show that w-lim n ^ oc T^ 1 / = w-lim n ^ oc Tfff. 
Indeed, let us choose r < min(t l5 Sj). From the inequality 

ll^^/llw = ll^w TLT nf\\n - \\ LT Hf\\n > 

which holds for t > r, we have 

||(Z#-2£)/|| w < k-tJ||L2?i/|| w -0, n^oo, 

which shows that T^™ / and T^™ / have the same weak limit. 

The set {T^/, t > 0} is bounded and hence sequentially weakly precompact, so that by 
passing to a subsequence, we may assume that 

w-limT*™/ = 7 , 7 e H . 

Next we show that for all (*) -sequences w-lim n ^ OG Tfy is a locally constant function 
fi — a.e. We begin by showing that 7 does not depend on the variables r. 

Let <9* ( „) denote the adjoint of d („) in 7i and let ^> G C^°(X). By the smoothness 

r i,m i,m 

properties of the density of the invariant measure we see that the function ip is in the domain of 
d* (v) and we have, using (3.27) 

(7,0V>VOw = lim (2£/,0; ( „,VOw = (d rW 7%f,il>) H = 0. 
Written explicitly, 

dpdqdrj(p,q,r)d M (il;(p,q,r)h(p,q,r)) = 0, (3.28) 

J i ' m 

for any ip E C£°(X). Since 7 e 7Y, we may set 7 = on the set A = {x E X ; h(x) = 0} and 
7 is locally integrable and thus defines a distribution in V'(X). By Eq.(3.28) the support of the 
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distribution d ( y ) 7(7?, q, r) does not intersect the set A and thus 7(7?, r) is fi-a.e.. independent 

r i,m 

of r. 

Let i > 0. Then w-lim n ^ 00 T^ +t f = T^y. Since t + t n | 00, it is easy to show, see 
[Br], that t n + t has an almost (*)-subsequence s n and from the above arguments we conclude 
that T^7 is independent of r. 

Next we show inductively, using Condition H2 that 7 does not depend on the variables p, 
q. Let (T^)* denote the semi-group dual to on H and denote Z its generator. 

Note that for ip e C^°(X) we have, upon integrating by parts 



= /dpdqdr £ [V rim • (V ri m + ftWi >m ) ^] 2* / 

m6{l,...,M} 

- J dpdqdr (L (^h))Ty, 



(3.29) 



where L is given by 

M M 

r „ . \7 -I- V n . ■ \7 - r\7 T/l • Y 

Pi 



m=l m=l j = l 



Since C^°(X) is in the domain of Z, we get 

A(^ 7 ) H = {Z^T*fi) H = limjZ^LT% +t f) n 
= (Z^,T^) n = - /dpdqdr (L (#))T^ 7 . 

The last equality follows from (3.29) since T.^7 is independent of r. 

We next choose ip(p, q, r) e C^°(X) of the form ip(r, p, 5) = <Pi(r)(p 2 (p, q)h~ l (p, q,r) 
with supp(<^ 1 (r)</? 2 (p, <z))nA = 0and Jdrip l (r) =0. For this choice of ip we have 

(^,T^ 7 ) W = Jdpdq(T^)(p,q)lp 2 (p,q)- JdrTp^r) = 0, 

and therefore 

= J dpdqdr 7 (p,q)L (7p 1 (r)lp 2 (p,q)) 

f f M 

= / dpdq 7 (p, q)V p ^ 2 {p, q)- J drJ2 r L, m Pi( r ) 

m=l 

f f M 

+ / dpdq 7 (p, q)V p Jp 2 {p, q) ■ J dr £ r R>m ^(r) . 

m=l 
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Since (p i (r) is arbitrary, it follows that 

J dpdq^{p,q)V p Tp 2 {p,q) = Jdpdq-f(p, q)V p Jp 2 {p, q) = 0, 

and thus, by a similar argument as above, ^(p, q) must be ^-a.e. independent of p { and p n : Thus 
7isafunction7(p 2 ,...,p n _ 1 ). 

Using this information, we choose now 

^(p,q,r) = (p l (p l ,p n )ip 2 (p 2 , ...,p n _ l ,q)^ 3 (r)h~ l (p,q,r) , 

with supp(ip Y ip 2 ip 3 ) fl A = and J dpjdp n ip l (p Y ,p n ) = 0. For such a choice of ip we obtain 

= J dpdqdrL (^ 1 ^ 3 )7 

= Jdp 2 ---dp n _ l dq 1 {V q ^ 2 )- Jdp l dp n p l ^ l Jdr7p 3 

+ Jdp 2 - • •dp n _ 1 dq 7 (V (?ji ^ 2 ) • J dp!dp n p n ^! J dxTp 3 , 

and from the arbitrariness of (p v (p 2 , (p 3 we conclude that 7 is independent of q l , q n (all our 
statements hold /j-a.e.). Finally, choose 

V>(p, q, r) = <p x (q l , q n )<fi 2 (P2, ■ ■ ■ >Pn-i » ?2» • • • » ?n-i WPi > Pn, (p, 9, r) , 

with supp^^y^) n A = and J dqjdq n ^(g^ g n ) = 0. Then we obtain 

= Jdp 2 . . . dp n _ 1 dq 2 . . . dq n-1 7(V p2 ^ 2 ) • J "dq 1 dq n (V q V)Tp x j d Pl dp n dr^ 3 

+ Jdp 2 . . •dp n _ 1 dq 2 . . .dq n _, l{V Pn _Tp 2 ) ■ J d qi dq n {V Qn _ 1 V)cp 1 J d Pl dp n dryJ 3 . 

From the arbitrariness of the (p i we conclude in particular that 

= |dp 2 . . .dp^dq, . . .dq n _ 1 7 (V p2 ^ 2 ) • J d qi dq n (V, V)^ . (3.30) 

We may choose (p l (q 1 ,q n ) = d („')<p(qi, q n ) for some v' e {1, . . . , d} and a positive (f>(q l7 q n ). 
By Condition H2 we see that 

X uy {q 2 ) = f d <li d <ln (d^V^^qJ = - J dq x dq n (d^d^Vfafa, q n ) 
is uniformly positive or negative. We can rewrite (3.30) as 

0= d P2--- d Pn-i dc i2--- dc in-i7d p MX Uy (q 2 )^ 2 , 

ve{\,...,d} 
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and we conclude that 7 is independent of p 2 . A similar argument shows that 7 is independent 
of q n _ 1 and iterating the above procedure we conclude that 7 is locally constant ^-a.e. 

So far, we have shown that for all (*)-sequences {t n } one has w-lim n ^ OG f = 7 = 
const. From the invariance of the measure /j, and its ergodicity we conclude that 

l=(hf)n = J»(dx)f(x). 

We conclude as in [Tr]: suppose that w-lim t ^ oc T^f 7^ (1, f) n . Then by the weak sequential 
precompactness of {T^f ; t > 0}, there exists a sequence u n | 00 for which w-lim^^ T^f = 
7] 7^ (1, f) n . But, referring again to [Br], the sequence {u n } has an almost (*) -subsequence 
{s n }. This implies that there is a (*)-sequence {t n } such that w-lim n ^ oc Tfyf = rj. This is a 
contradiction, since we have seen that w-lim n ^ OG Tfy-f = (1, for all (*)-sequences. By a 
simple density argument this implies that 

t i™ / / u ( dx )/( x ) T 'Hf <» = y /"( dx )/( a; ) y /«(dx)^(x) , 

for all f,g ETC and the proof of Proposition 3.9 is complete. 

With Proposition 3.9 the proof of Theorem 2.2 is now complete. 



4. Commutator Estimates and Spectral Properties of 

In this section, we prove Proposition 3.4 and Proposition 3.5. We generalize the commu- 
tator method of Hormander to study the spectral properties of the operator L Hq which is, by 
Lemma 3.1, the closure of the differential operator L with domain C£°(X) which we defined in 
Eq.(3.11). We recall the definition: 

M >2 



m=l 
M x 2 



+ E ^T^^ V PRWRm y V (4.1) 

m=l ^ V 7 

M M 



with the abbreviations 

W L,m = A L,m r L,m ~ 9l , W R,m = K^^m ~ Qn , ( 4 - 2 ) 

and where L s is the Liouville operator associated with the Hamiltonian H s (q,p): 

n 

L s = E^-V 9j -(V g ^).V pj . (4.3) 

3 = 1 
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For the following estimates it will be convenient to work in the flat Hilbert space L 2 (X, dx) . 
The differential operator L is unitarily equivalent to the operator /C on L 2 (X, dx) with domain 
C °°(X) given by 



K = e -PoG/2 Le P G/2 

M \2 M \ 2 (A A) 
V" L,ro 'L,m p * p \^ A R,m7R,m pH< r l^'^ 

m=l Pl m=l P R 



where a is given by (3.10) and 



= V rL|m + J (A -P /2)P /2W Um , 



«R,m = V rR m + yJilk-Po/lWlW^ , 

M M 

= 1 







M 

E 

m=l 


A 2 T 

L,m 'L,m 


& 


+ /3r 


2 


& 


-/3r 


M 

E 

m=l 


A R,m7R,m 


& 


+ Pr 


2 



All subsequent estimates will be valid for any / e 5(X) and thus for all functions in the domain 
of/C. 

It is convenient to introduce the following notations: We introduce new variables, and 
recall some earlier definitions: Let n = [n/2] denote the integer part of n/2. We define 





= V p . + a LPj , j 


= l,...,n' 


i 




Pj 


= + a R pj , j 


= n' + l,.. 


.,n , 




Qi 


= ^ qj +a h W J (q,r), 


.7 = 1,. 


. . , n . 




Qi 


= ^ qj +a R W j (q,r), 


J = w' 


+ 1,.. 


.,71 , 




= V r L , m + a L W L,m(Q> 


r) , m 


= 1,.. 


.,M, 


-^R,m 


= V r R , m + a R W R,m(^ 


r) , m 


= 1,.. 


.,M, 



Anharmonic Chains Coupled to Two Heat Baths 



25 



where 



a T = 



((/3 L -/3 /2)/3 /2) 1/2 
% = ((/3 R -A)/2)/y2) 1/2 



M 



L,m ' 



m=l 



W n (q,r) = V q .V(q), j = 2,...,n-l, 



M 



W n (q,r) = V qn V(q)-J2 r R,m, 



W Um(<l, r ) = X L 2 m r L,m ~ ?1 . TO = 1, . . . , M , 



V(«' r ) = Vm r R,m - 9 n , TO = 1, . . . , M . 

We next define the operators K , K, and A which will be used in the statement of our 
main bound: 



1 + E P I P j + E QjQj + E ( R t,m R L,m + R R,m R R,m) ) 
j = l j = l m=l ' 

M 

K = a — K, = -Kq + ^ ( b L,m R L,m R L,m + b R,m R R,m R R,m) • 

Here, we use 



1/2 



(4.5) 



m=l 



J5 m 

Our main estimate is 



b h,m _ \,m7L,m/ /^L ) b R,m ~ \,mlR,m/ Pr • 



Theorem 4.1. Under the Assumptions HI, H2 on V, there are an £ > and aC < oo such 
that for aJJ / e one has 



||A e /|| < c(||k/|| + n/ii). 

Proof. The proof will be an easy consequence of the following 



(4.6) 



Proposition 4.2. There are finite constants C-, C'-, and C such that for all f e S(X) one 
has with n = [n/2], 



\\A^- l PJ\ 



W^-'Pn+l-jfl 



\\^- l Qn +J -J\ 
\ R L,mf\\ + \\ R R,mf\ 



< C,(||K/|| + ll/H), J = l,. 

< + 11/11), 3 = h- 

< C j ()\Kf\\ + 11/11) , J = h- 



< c;(iik/h + n/ii), j = i,. 

< C(\\Kf\\ + H/ll), to=1,. 



,ri , 



. , n — n , 
. ,n' , 
. ,n — n' , 
.,M, 



(4.7) 
(4.8) 
(4.9) 
(4.10) 
(4.11) 
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where e- = 4 1 2j and e'- = 4 2j . 

Proof of Proposition 4.2. For the R { m , we have the easy estimate 

\\R hm f\\ 2 = (f,KmRi, m f) < KXx*(f,Kf) < b-X\\Kf\\ ii/n < b-liWKfW + Wfwf. 

(4.12) 

This proves Eq.(4.1 1) for these cases. 

For the other cases, the proof will proceed by induction: It will proceed by bounds on P 1 , 
Q l , P 2 , ■ ■ ■ , Q n ,, and a totally symmetric argument, which is left to the reader, can be used from 
the other end of the chain, proceeding over P n , Q n , P n _ l , until the bounds reach the "center" 
of the chain. We next prepare the inductive proof. To make the result of this calculation clearer, 
we define the matrices 

M 3,k = V q 3 V i k V > j,k = l,...,n. 
In components, this means, for /x, v e {1, . . . , d}, 

M jT ] = V„wV>)7, j,k=l,...,n. 
By our choice of potential V all the Mj k vanish, except M.- ., withj = 1, . . . , n and -M J+1 ■ = 

■Mj with j = 1, . . . , n — 1. Furthermore, by Condition HI, all the M. ?^ are uniformly 
bounded functions of q. Finally, by Assumption H2, the matrices M.. . +1 are definite, with 
uniformly bounded inverse. 

One verifies easily the relations: 

[P^K ] = Q,, 

M 

[Q V K ] = -M l l P l - M 2 l P 2 + ^2 C 'm( R L,m + R L,m) > (4-13) 

m=l 

[P p K Q ] = Qj , j = 2,...,n-l, 

[g,,K ] = M, ,.,/> , M^I) M j+hj P j+l , , 2 n 1. 

where 

C m = 7 L ,m(/?L - Pq)/Pl > 
C m = ^^L^L ~ A)) • 

Symmetrical relations hold at the other end of the chain. With these notations, we can rewrite 
(among several possibilities): 

Qy = [P 1; K ], 

M 

p 2 = -M 2 -;([Q„i( ]+Vi-E c ™^+U' ( 4 - 14 ) 

m=l 

Q J = [P,,^], j = 2,...,n, 
P j+ i = -MjlrjdQ^KJ + M^P^ + M^) , j = 2,...,n', 
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with symmetrical relations at the other end of the chain. We can streamline this representation 
by defining Q = R Ll , and M 10 = — 1. Then we can write, for j = 1, . . . , n: 

Pj = -Mj^ilQ^KJ + Sj) , (4.15) 
Qj = [P^Kq], (4.16) 

where the operators <S - depend linearly on {P 1: . . . , -P,_i}, {<5i, • • • , Qj_\}, and the R L m . The 
relations Eqs.(4.15) and (4.16) will be used in the inductive proof. 

Such relations are of course reminiscent of those appearing in the study of hypoelliptic 
operators. The novelty here will be that we obtain bounds which are valid not only in a compact 
domain, but in the unbounded domain of the p's and q's. 

The following bounds will be used repeatedly: 

Proposition 4.3. Let Z denote one of the operators Qj, Q*, Pj, or P* . Let M. denote one 
of the Mj k . Assume that a G (0, 2) . Then the following operators are bounded in L 2 (X, dx) : 

1) A ,/3 LM,A~ a ]A 7 , if P + l < a+ 1 , 

2) t^ZK 1 , if/3 + 7 < -1 , 

3) A^K^ZW ,ifp + j<-l, 

4) A^A-^A 7 ,if/3 + 7 < a+ 1 , 

5) A^A-^i^A 7 , if 13 + 7 <a. 

Proof. The proof will be given in Appendix B. 

Because we are working in an infinite domain, and work with non-linear couplings, we 
will not bound the l.h.s of Eq.(4.7) directly, but instead the more convenient quantity 2 : 

Rj(f) = ia-' ,/y . A- '/>/). 

We have the 

Lemma 4.4. There is a constant C such that for all j E { 1 , . . . , n} and all f G S(X) one has 
the inequality 

W'-^fW 2 < cdR^i + Wff). 

Therefore, to prove Eq.(4.7), it suffices to prove the corresponding inequality for the Rj. 

Proof of Lemma 4.4. Let M = Mj -_ l5 e = e -, and P = Pj. Then, by our Assumption 
H2, there is a constant m > for which M. > m. Therefore, 

IIA^iVH 2 = (A e ~ l Pf , A e ~ l Pf) 

< m- l \{MA e - l Pf , A e_1 P/)| 

< m- l \(A e - l MPf , A £ - l Pf)\+m- l \([A £ -\M]Pf, A e ~ l Pf)\ 

< m- 1 | J R J (/)| + m- 1 |((A £ [A £ - 1 ,A^]A)(A- 1 P)/, (A~ l P)f)\ . 

2 

For readers familiar with the method of Hormander, we wish to point out that this device seemed necessary 
because we do not have good bounds on [K Q , [Q l ,K Q ]]. 
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The proof of Lemma 4.4 is completed by using the bounds 1) and 2) of Proposition 4.3. 

The inductive step. We begin by the induction step for the P-. We assume now that the 
bounds (4.7) and (4.9) have been shown for all j < k. We want to show (4.7) for j = k + 1. 
Using Eq.(4.15) and Lemma 4.4, we start by writing 

R k+ i(f) = {^ ek+l -'M k+hk P k+ J , A^W) 
= ^A 2ek+1 ~ l [K , Q k ]f , A~ l P k+l f^j 

= X l - x 2 . 

We first bound X 2 . Note that S k+i is a sum of terms of the form MT where T is equal to P- or 
with j < k, and A4 is either a constant or equal to one of the M k e . Therefore we obtain, 
using Proposition 4.3, the inductive hypothesis, and the choice 2e k+l < min - <fc (e ■, e'-) = e' k : 

\(A 2e ^- l MTf , A~ l P k+ J)\ 

< \{MA 2e ^~ l Tf, A- l P k+ j)\ + \{([A 2 ^-\M]A)(A- l T)f , A~ l P k+l f)\ 

< 0(l)(||iv-/|| + ||/||)||/||+0(l)||/|| 2 < 0(l)(||iv-/|| + ||/||) 2 . 

This proves the desired bound. 

We now come to the "interesting" term X v The commutator is rewritten as 

[K , Q k ] = -Q k K - K*Q k + \ (Q k (K + K*) + (K + K*)Q k ) 

= x 3 + x 4 + x 5 . 

We discuss the 3 corresponding bounds: 

Term X 3 . In this case, we are led to bound, with e = e k+1 , 

T 3 = \(Q k Kf, A 2s ~ 2 P k+ J)\ = \(Kf , QlA 2s - 2 P k+ J)\ 
= \(Kf, (Q* k A 2s ~ 1 )(A~ 1 P k+1 )f)\ 

< \(Kf, (A-^XQJA 2 *- 1 )/)! + VQl^ 2e -\^ l P k+ m 

= X 3,l+ X 3,2 ■ 



(4.17) 



We start by bounding X 3 { . Since A 1 P k+l is bounded by Proposition 4.3, it suffices to show 
that 

WQt^-'fW < C(\\Kf\\ + H/ll). (4.18) 
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To see this we first write, using Q = Q k , 

||g*A 2e - 1 /H 2 = (/ , A 2£ - l QQ*A l£ - l f) 

= WA^Qff+if, [A 2e ~ l Q, Q*A 2e ~ 1 ]f) . 

The first term is bounded by the inductive hypothesis by 0(1) + ||/||) and the choice of 

e k+l , while the second can be bounded by 0(1)||/|| 2 by expanding the commutator (and using 
Proposition 4.3): 

[A 2£ - 1 Q,Q*A 2£ - 1 ] = (A 2£ - 1 Q*A- 2e )A 2e [Q,A 2e - 1 ] 
+A 2e " 1 [Q, Q*]A 2e_1 + ([A 2s ~\ Q*]A 2e ) A~ l Q . 

This proves Eq.(4.18). 

To bound X 3 2 , we use [P£ +l , Q k ] = and we write 

[QlA 2s ~\A-'P k+l ] = (<5fcA _1 )[A 2e_1 , P fc+1 ] + ([Q^, A _1 ]A 2e ) (A~ 2e P fc+1 A 2e_1 ) . 
Since each factor above is bounded by Proposition 4.3, the desired bound follows: 

>2 



T 3 < o(i)(||*:/|| + 



Term X 4 . Here, we want to bound T 4 = \{K*Q k f , A Z£ - z P k+ J)\. We get 

T 4 = \{K*Q k f, A 2s ~ 2 P k+ J)\ = \(Q k f , KA 2e ~ 2 P k+l f)\ 
< \(A 2e ~ l Q k f , A~ l P k+l Kf)\ + \(Q k f , [K,A 2e " 2 P fc+1 ]/)| (4.19) 

= X 4 l + X 4j2 . 

Using the inductive hypothesis, and the bound ]| A _1 P fc+1 1| < 0(1), the term X 4 j is bounded 

by 

W^-'QJW W^P^KfW < o(i)(||K/|| + ||/||) 2 . 

We write the commutator of X A 2 as 

[K,A 2s ~ 2 P k+l ] = A 2s - l (A- l [K,P k+l ] + A- l [K ,A 2 - 2s ]A 2s - l (A- l P k+l )) , 

since K — K commutes with A. Using Proposition 4.3 and the inductive hypothesis this leads 
to the following bound for X 4 2 : 

X 4a < \(A 2e ~ l Q k f , A~ l [K, P k+l ]f)\ 

+ \ {^ 2e - l Q k f, (A- 1 [K ,A 2 - 2e ]A 2e - 1 )(A- 1 P fc+1 )/)| 
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This completes the bounds involving X 4 . 
Term X 5 . Here, we bound 

T 5 = i((Q k (K + K*) + (K + K*)Q k )f, A 2e " 2 P fc+1 /) . 

Assume first k > 1 (and in any case we have k < n). Looking at the definition of K, we see 
that in this case Q k commutes with j(K + K*) = Re K, and we can rewrite T 5 as 

T's = 2^^/, Q*A 2e " 2 P fc+1 /) . 

Using the Schwarz inequality and the positivity of Re K , we get a bound 

|T 5 '| < ((ReK)/, /) 1/2 ((ReK)Q*A 2e - 2 P, +1 /, g*A 2£ " 2 P fc+1 /) 1/2 
= (Rc(Kf , f)) l/ \Rc(KQlA 2s - 2 P k+ J , QlA 2s - 2 P k+ J)Y /2 
= (Re(Kf , f)) l/2 (Re(A~ 2s KQ* k A 2e ~ 2 P k+l f , A 2s Q* k A 2s ~ 2 P k+l f)) 1/2 
ee (Re (Kf , /))^ 2 (Re (f l , f 2 )) 1 ^ 2 • 

1 /2 

The first factor is clearly bounded by (||X/|| + || / 1| ) .To bound f l , we expand again: 

f t = A~ 2e KQ*A 2£ - 2 Pf = (A- 2e Q*A 2e - l )(A- l P)Kf 
+ A- 2e [K,Q*]A 2£ ~ 2 Pf + A~ 2e Q*[K,A 2£ - 2 ]Pf 
+ A- 2e Q*A 2£ ~ 2 [K,P]f . 

The norm of the first term is bounded by 0(1) + ||/||). Using Proposition 4.3, the other 

terms are bounded by 0(1)||/||. To bound f 2 we write 

f 2 = A 2e Q*A 2e ~ 2 Pf = A" 1 PQ*A 4£ - 1 / 
+ A 2e Q*A- 2£ - 1 [A 4e - 1 ,P]/ 
+ A 2e [Q*,A- 2e - 1 ]PA 4e - 1 / . 

We control the first term using the inductive hypothesis (it is here that we use the factor 
4e < e' k ) and the two others by Proposition 4.3. Combining these bounds, we finally get the 

bound T 5 < 0(1)(||K/|| + ||/||) 2 , and hence the inequality (4.7) is shown for all j. 

It remains to discuss the cases k = 0, 1 for the term X 5 . The commutators of Re K with 
Q = i? L j or with Qj do not vanish and hence there are additional terms in T' 5 . They are of the 
form 

M 

^2 b L,m([ R L,m R L,mi R L,l^f » A 2e_2 P fc+1 /) , 
m=l 

M 

E KmHRtmR^ Qllf , A 2£ " 2 P fc+1 /) . 
m=l 
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Since [#L,rA,m> #L,l] = COnst -#L,l 5 m,l and \-K,m R L,m, Ql\ = COnst -#L,m> this is obvi " 

ously bounded by 0(1)(||K/|| + 



We have discussed now all the cases for the inductive bound on the P . The discussion of 
this step for the Qj is the same, except that some simplifications appear because of the simpler 
relations Q- = [P J7 K ]. The proof of Proposition 4.2 is complete. 

Proof of Theorem 4.1. Lete < e' n+1 . We rewrite 

n+l n 

A 2e = A 2e " 2 (i + Y t Q* j Q j + J2 p j p j) ■ ( 4 - 2 °) 

3=0 3=1 



Note now that for Q = Q., 

A 2e ~ 2 Q*Q = Q*A 2e ~ 2 Q + [A 2s ~\Q \Q 
Using Proposition 4.2 and Proposition 4.3, we get a bound 

(f , A 2e ~ 2 Q*Qf) < 0(i)(||iv-/|| + ||/||) 2 + 0(i) 

Of course, the P satisfy analogous relations. Since ||A 2e_2 /|| < 0(1)||/||, the assertion (4.6) 
follows by summing the terms in Eq.(4.20). The proof of Theorem 4.1 is complete. 
Using Theorem 4.1 we can now prove Proposition 3.4. We have 

Proposition 4.5. If the potential V satisfies Conditions HI, H2 and if (3 < 2min(,9 L , /3 R ) 
there is a A* > such that if the couplings satisfy |A L |, |A Rm | G (0, A*) then both L Ho and 
L* H() have compact resolvent. 

Proof. We show that the operator /C on L 2 (X, dx) has compact resolvent. From Theorem 4.1 
we get the bound 

||A e /ll < C(||(/C-a-l)/|| + ||/||), (4.21) 

for all / G S(X). Since, by Lemma 3.1, C^°(X) is a core of /C, we see, by taking limits, that 
the estimate (4.21) holds for all / in D(K). 

We note that A 2 has compact resolvent. Indeed, recall the definition Eq.(3.7) of the 
effective Hamiltonian G. It is easily checked that, first of all, G grows quadratically in every 
direction of R 2,i ( n + M ) 5 f or sufficiently small |A i>m |. Second, it is also easily verified that 

n M 

A2 = 1 - B A W + V + E + A r R ,J + W(p, q, r) , 

j = l m=l 

and, by construction 

n' n 

W(p, ? ,r) ^ ^a 2 ((V pj G) 2 + (V Qj G) 2 ) + ^ a 2 ((V pj G) 2 + (V gj G) 2 ) 

j=l j=n' + l 

M 

+ E(«L(V ri ,„G) 2 + 4(V rR ,„G) 2 ) , 



m=l 
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up to bounded terms. Thus W(p, q, r) diverges in all directions of R ' . Using the Rellich 

criterion (see [RS], Thm XII.67) we conclude that A e has compact resolvent for every e > 0. 

Therefore, Eq.(4.21) implies, using again the Rellich criterion, that (JC — a — 1)*(/C — 
a — 1) has compact resolvent. We claim this implies that JC itself has compact resolvent. 
Indeed, since JC — a — 1 is strictly m-accretive, its inverse exists, and therefore the operator 

((JC-a- \)*{JC-a- l))" 1 = ((JC-a- l)*) _1 (/C-«- l)" 1 exists and is compact. This 
implies that (JC — a — 1) _1 is compact and hence JC has compact resolvent as asserted. 
Finally, we prove Proposition 3.5. We have the following 

Proposition 4.6. Let g denote an eigenvector of L Ha or . If the assumptions of 
Proposition 3.4 are satisfied then gexp(/3 G/2) is in the Schwartz space S(X). 

Proof. We prove the corresponding statement for the operator JC on L 2 (A, dx). We consider 
the set of C°° vectors of e Kt , i.e the set 

C°°(JC) = {/ G L 2 (A,dx) ; e Kt G C°°(R + , L 2 (A, dx))} . 

The set C°°(JC) obviously contains all eigenvectors of JC. Therefore Proposition 4.6 is a direct 
consequence of the following Proposition. 

Proposition 4.7. C°°(JC) = S(X). 

Proof. By Theorem 1.43 in [Da] we have the following characterization of C°°(JC): 

C°°(fC) = n n > D(JC n ) , 

where D(JC n ) = {/ € D(JC n - 1 ), JC n ~ l f G D(JC)}. 

Since S(X) C D(JC) and JCS(X) C S(X), we have the easy inclusion 

S(X) c n n D(JC n ) = C°°(JC) . 

To show the inclusion in the other direction we will need the following Theorem which we 
will prove in Appendix C. This is a (slight) generalization of the core Theorem, [Da], Thm 1.9. 

Theorem 4.8. Let B be a Banach space. Let A : D(A) — > B be m-accretive. For all 
n = 1 , 2, . . ., if V is a subset of D(A n ) and is dense in B and furthermore V is invariant under 
the semi-group e At , then V is a core for A n . 

Given this result we first show that S(X) is invariant under e *. For s > we consider 
the scale of spaces J\f s given by J\f s = D(A S ), with the norm ||/||( s ) = ||A S /||. Fors <0welet 

N s be the dual of M_ a . From the definition of A , it is easy to see that {|| • ||( s ) ; s = 0, 1, • • •} 
is a system of semi-norms for the topology of S(X) and hence S(X) = n s A/" s . 

To show that S(X) is left invariant by the semi-group e Kt generated by JC, it is enough to 
show that 

e Kt N s C N 8 for all s > . 
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For /, g in S(X) we have the identity 

(A-e^A'/ ,g) = {f, A'e&A-'g) 

= (f,g)+ [ dr(/, A s JCe^A- s g) 
Jo 

= (f,g)+ f dr{f, (K, + B)A 8 e 1CT K~ 8 g) 
Jo 

= (f,g)+ [ dr(A" s e K;V A s (/C* + B*)f , g) , 
Jo 

where 

B = [A s , /C]A" S , 



(4.22) 



is a bounded operator by Proposition 4.3. From (4.22) we see that 

-^-A^e^A 5 / = A-'e^A't/C* + B*)f . (4.23) 



Now /C* is the generator of a strongly continuous quasi-bounded semi-group, B* is bounded and 
so, [Ka], Chap. 9, Thm 2.7, /C* + B* with domain D(IC*) is the generator a strongly continuous 
quasi-bounded semi-group e (/c * +ir) * with ||e (/c * +B * )t || < e (a+l|s * From (4.23) we see 
that 

e (K*+B*)t = ^s e !C*t A s _ 

Thus we obtain 

|| A - e /C-t A *|| < e («+||B'||)t ; 

and so e' ** : A"_ s -> A"_ s , s > 0, is bounded. By duality e Kt : N s ->• A" s , s > 0, is also 
bounded. This implies that 

a; c N 8 , s > , 

and therefore 5(A) is invariant under e Kt . 

We now use Theorem 4.1. Let / e 5(A), then replacing / by A m / in Eq.(4.6), we obtain 

ll/IWe) < 0(1) (l|/CA m /|| + ||/|| (m) ) 

< o(i)(||/c/|| (m) + ||[/ir,A m ]/|| + ||/|| (m) ) . 

Since 

||[/C,A™]/|| = ||A-[/C,A— ]A m /||, 
and since A m [/C, A~ m ] is bounded by Proposition 4.3 we obtain the bound 

(m+e) < 0(l)(||£/|| (m) + ||/|| (m) ) . (4.24) 
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Using (4.24) it is easy to see, by induction, that, for n = 1, 2, • • • we have 



(ne) 

Since is a core for /C n by Theorem 4.8, we see, by taking limits, that 

D{K?) C Af n£ ■ 

Therefore 

C°°(JC) = n n D(iC n ) c n n Af ne = S(X) . 

And this concludes the proof of Proposition 4.7. 
Appendix A: Proof of Lemma 3.1 

If x(t) = £(£, w; x) denotes the solution of (3.1), it has the cocycle property 

£(£, t s w; f (s, w;x)) = £(t + s, w; x) , 

which holds for all t, s £ R, x G X and w G W. Here we have introduced the shift 
(r t it;)(s) = w(t + s) on W. In particular the map x h- > is a bijection with inverse 

x i— > £(— £, r^; x). A standard argument shows that these maps are actually diffeomorphisms 
(see e.g., [IW], Ch. V.2). The Jacobian of w; ■) is given by 

J(t,w;-) = \detD x t(t,w;-)\ = '^"^ , 

and according to (3.2) the Jacobian satisfies 

< J(t,w;x) < e B l*l . 

Remark. In our case we have in fact 

div6 = -dJ2% m = -r<o, 



so that 

J(t,w;-) = e~ Tt . 
The Lemma 3.1 is an immediate consequence of the following lemmata. 
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Lemma A.l. T extends to a strongly continuous, quasi-bounded semi-group of positivity 
preserving operators on L 2 (X, dx). Its generator is the closure of L. 

Proof. Let / eC °°, then we have 

IITVII 2 = lim f dx|TV(x)| 2 = lim [dx X{]x{<R} |E[/ o £(£, w; x)} | 2 

" jft^o / dX ^{|x|<i?} E 0/| 2 °^ t '« ; ^)] 
- ^ E [J d *X{\x\<R}\f\ 2 °Z( t > W ' X )] 

^ R li Zo E t/ dy *{M<*} ° y) l/|2(?/) J(_t ' y) ] 

< lim E[[dyl\f\\y)e m ] = e m \\ff . 

Thus T* extends to a bounded operator on L 2 (X, dx) by continuity. A simple approximation 
argument shows that T* is weakly continuous at £ = 0, and hence, since it is obviously a 
semi-group, strongly continuous. Positivity is evident. Now let / e C™(X), then 

[(l/f)(T* - 1)/ - L/](x) = T d,(l/£)(T s - l)Lf(x) , 

Jo 

from which we conclude that the generator L of T* on L 2 (X, dx) satisfies L C L. 

From the inequality Re (/, L/) = -±\\a T \7f\\ 2 - (/ , div bf) < B\\f\\ 2 , and the fact 
thatC °° C D(L*) and Re (/, L*/) = -\\\a T \7f\\ 2 + (/ , div bf) < B\\f\\ 2 , one concludes 
that: 

(i) L is accretive. 

(ii) The range of (A — L) is dense for Re (A) > B. 

Hence, by the Lumer-Phillips Theorem (see for example [Da], Theorem 2.25), the closure 
L generates a quasi-bounded semi-group on L (X, dx). Since such generators are maximal 
accretive, we conclude that L = L. 

We shall now consider the Markov semi-group on weighted L 2 -spaces (such as H Q ) of the 

form 

H s = L 2 (X,e" s dx) , 
where S G C°°(X), and e~ s eL^dx) is normalized (He" 5 ^ = 1). We also assume that 

b s = DVS 

with D as in Eq.(3.3) satisfies the condition 

Hdiv&sIL < oo . (A.l) 
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The action of T* on H s is obviously equivalent to that of e~ s,2 T t e s/1 on L 2 {X, dx) . For 
/ G C °°(X), Ito's formula gives 

{e S/2 T t e S/2 )f{x) = E [ e (S(Ht))-S(H0)))/2 f{xmx{0)= ^ 

= E[e^r ds ( iS )(^ s )) + ^o td -( s )^ Tvs )(^ s ))/ (a; (t))|x(0) = x] 
= E[D(t)e^o dsRs(x(s)) f(x(t))\x(0) = x] , 

where 

R s (x) = (LS)(x) + ^(VS-DVS)(x) = div (b s ) + (b + l -b s ) ■ VS . 
By the Girsanov formula we obtain 

{e -S/2 T t e S/2 f){x) = E[e U> R s(v(s)) f{ymm = x] 

where y(t) is the Markovian diffusion process defined by the equation 

y(t) = y(0)+ [ ds(b + b s )(y(s))+a(w(t)-w(0)). 
Jo 

Assuming that R s is bounded above: 

E s = sup R s (x) < oo , 
xex 

and denoting by T l s the semi-group on L 2 (X, dx) associated with the process y(t), by Lemma 
A.l we get 

\ {e S/2 T t e S/2 f){x)l < e ^st {T t mxh 

from which one concludes that e~ s ^ 2 T t e s ^ 2 extends to a strongly continuous, quasi-bounded 
semi-group of positivity preserving operators on L (X, dx). 

By the Feynman-Kac formula (or Cameron-Martin) we can conclude that the generator of 
this semi-group is given, on C™(X), by 

L s = V • DV + (b + b s ) ■ V + l -R s = e~ s/2 Le s/2 . 

Repeating the argument of Lemma A.l we conclude that L s is the generator. Since C^°(X) is 
invariant by e ±s ^ 2 we get 

Lemma A. 2. Let S £ C°°(X) be such that 
(i) b s = DVS satisfies Condition (A.l), 
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(ii) sup xex (6 + \b s ) ■ VS(x) < oo. 
Then the semi-group T* extends to a strongly continuous quasi-bounded semi-group on H s . 
Moreover (X) is a core for its generator. 

Now Lemma 3. 1 is a direct consequence of Lemma A.2. Indeed we apply Lemma A.2 to 

S(x) = p G(p,q,r), 

where G is given by (3.7). We see that Condition (i) of Lemma A.2 is satisfied, since G(p, q, r) 
is of the form quadratic + bounded. An explicit computation shows that the assumption 

/3 < 2min(/3 L ,/3 R ) , 

implies that Condition (ii) of Lemma A.2 is satisfied and that the semi-group satisfies the bound 
||T*|| < e at where a is given by Eq.(3.10). 

Appendix B: Proof of Proposition 4.3 

To prove the claims it is useful to introduce some machinery which replaces the pseudo- 
differential calculus, which seems unavailable for the class of operators we want to consider. 
This may be useful in its own right. 

Let J 7 , as in the Hypotheses HI, H2 denote the class of functions of q G R dn which are 
bounded together with all their derivatives. Let y denote the linear space of operators spanned 
by 

f(q)q m dfp n d;'r e d £ r ' , (B.l) 

where / G T . (The notation is sloppy, we really mean components v = 1 , . . . , d of each of these 
quantities.) We shall say that the quantities in Eq.(B. 1) are of degree m + m + n + n + £ + £'. 
We let y s denote the subspace of y spanned by the expressions of degree less than or equal to 
s. Clearly, the operators Z and M of Proposition 4.3 are in y l and y°, respectively, while K Q 
and A are in y . Below, we shall use this, but also an additional property of the potential V. 
We have 

Lemma B.l. One has the inclusion 

[A 2 ,y s ] c y s+l . (b.2) 

Furthermore, Y G y° defines a bounded operator. 
Proof. By inspection. 

Proposition B.2. Assume thatY G y J , for some j G {0, 1, . . .}. Then 

A^yA 7 (B.3) 
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defines a bounded operator on L 2 (X, dx), when 

P + 7 < -3 ■ 

Let Z be an operator in y. Assume that [A 2 , Z] G y\ for some j G {0, 1, . . .}. Then 

A^[A- a ,ZW (B.4) 
defines a bounded operator on L 2 (X, dx) for all a, j3 and 7 satisfying 

/3 + 7 < a-j + 2. 

We will give bounds on various quantities involving y s . For this, we will use throughout 
the following device: 

Lemma B.3. Let A z be a bounded continuous operator-valued function of z and let F(\, z) 
be a real, positive continuous bounded function. Then 

poo poo 

II / dzA z F(A,z)u\\ < sup \\A \\ \\u\\ / d2supF(A,^). (B.5) 

Jo ' y>0 Jo A>1 

If furthermore A = A z is independent of z, one has the bound 

POO POO 

II / dzAF(A,z)u\\ < \\A\\ H^ll sup / dzF(\,z). (B.6) 
Jo \>\Jo 

Proof. Note first that 

POO POO 

11/ dzA z F(A,z)u\\ < / dz\\A z \\\\F(A,z)u\\ . 
Jo Jo 

Since A, as an operator, satisfies A > 1 we also have from the spectral theorem: 

||F(A,z)|| < supF(A^). 

A>1 

Thus, Eq.(B.5) follows. In a similar way: 

POO POO POO 

11/ dzAF(A,z)u\\ < \\A\\ || / dzF(A,z)\\ \\u\\ < \\A\\ \\u\\sup / dzF(X,z), 

Jo Jo A>1 JO 

which is (B.6). The proof of Lemma B.3 is complete. 

We shall also make use of the following identity, valid for a G (0, 2), [Ka] Thm. V.3: 

A _ Q= SWW2) r dzz -a/2 {z + A 2 r ^ {BJ) 
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We also let C_ a = sin(7ra/2)/7r. 

Proof of Proposition B.2. It is obvious that if we show the claim for (3 + 7 = —j, then it 
also follows for (3 + 7 < —j. By the definition of y 3 , and observing that f(q) is bounded, and 
by the explicit form of A , we see that the claim holds when 7 = 0. We next consider the case 
/3 + 7 < 0, — 1 < 7 < 0. In this case we write 

A /3 yA 7 = YA 13 ^ + A^[Y, A 7 ] . 

The first term is clearly bounded as in the case 7 = 0, by considering adjoints. The second term 
can be written as 



/•OO 

A?\Y, A 7 ] = C / Azz^h? 
Jo 

JO 



y 



1 



z + A 2 



A/3 

dz Z T/2_ii__ [ y > A 2 ]- 



1 



(B.8) 



, z + A 2L ' J ^ + A 2 ' 

By Lemma B.l, we see that [Y, A 2 ] e y 7-1-1 and thus, we get, using Eq.(B.5), 



||A^[Y,A 7 ]|| < C 7 sup| 

y>0 



A /3 



y + A 2 



/•oo 1 

[Y,A 2 ]|| / dzz^sup 

JO AM 2 + 



A 2 



The norm is bounded because 7 e [—1,0) and thus /3 — 2 + j + l < /3 + 7 + j = 0. The sup 
over A is (1 + z) - and the integral converges because 7 e [—1,0). 

We now proceed to the other choices of 7 by induction. We first deal with negative 
7. Assume we have shown that A^YA 7 is bounded for all 7' e [— r, 0], and assume that 
7 e [-r - 1, -r), and Y e y 3 . We write 

A^YA 7 = A^ _1 yA 7+1 + A /3 [Y,A~ 1 ]A 7+1 . 

The first term is bounded by the inductive hypothesis. To bound the second, we apply again the 
method used in Eq.(B.8). Then we get 



r°° a/3 
A^[Y, A _1 ]A 7+1 = C_, / dz z~ l/2 ——r[Y, A 2 ] 

Jo ^ + A 

Since [Y, A 2 ] e we see from the inductive hypothesis that 



A 7+1 
z + A 2 



sup 

y>0 



A 13 



y + A 2 



[Y, A 2 ]A 7+1 



< 



A P-2 [Yj A 2 ] A 7+1 



is bounded and hence, using (B.5), we can complete the inductive step. 

The case 7 > is handled by observing that Y E y 3 implies Y* E y J , and bounding 
A' 3 YA 7 by bounding A 1 Y* A^ . This completes the proof of the first part of Proposition B.2. 
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To prove the second part, we first consider the case a G (0,2). In this case, using Eq.(B. 7), 



we write A /3 [A" Q ,Z]A' 



C „ / d^^- Q/2 A^ 



as 



+ A 2 



'°° \@ AT 

\ - C_„ I dzz~^ 2 -^ [A 2 , Z] ' 



-a 







+ A 2 



z + A 2 ' 
(B.9) 



We let B = [A , Z], use another commutator and rewrite (B.9) as 



' A 7 /"°° \P A 7 

C_ a / dzz- a ' 2 hPB - t \ 2 , 2 + C_ a / d^-/ 2 ^— [A 2 ,/j1 ^ 

'„ (z + A 2 ) 2 a J z + A 21 j (z + A 2 ) 2 



= C_JX 1+ X 2 ). 
We first bound X v We get, using Eq.(B.6), 



, 1 A 3 

(z + A 2 ) 2 



< IN 

< IM 



A /3 BA -J-/3 

A^BA-i- 13 



yj+P+i 



5M Az ^' 2 VTW- 



sup 

A>1 70 



d SS -"/ 2 A 2 - Q 



yj+P+i 
(s+ 1) 2 A 4 



Since, by assumption, £? e y 7 , the norm is bounded, and the integral is bounded because 
P + 7 < a — j • + 2, by assumption. 

To bound X 2 we first observe that by assumption, and by Lemma B.l, C = [A 2 , B] is in 
y-* +1 . Therefore, using Eq.(B.5) we find the following bound for X 2 : 



\X 2 u\\ = 



o A/3 AJ-1+/3+7 



z + A 2 



(z + A 2 ) 2 



< \\u\\ sup 
y>o 



A 13 



y + A 2 



dzz~ Q / 2 sup 



o 



xVi (z + A 2 ) 2 



This is clearly bounded when /3 + 7<a — j + 2 and a £ (0, 2). This completes the second 
part of Proposition B.2 when a £ (0, 2). If a = the assertion is trivial. The case a = 2 is 
handled by considering the identity: 



[A~\Z] = A" 2 [A 2 ,Z]A _/ . 
The cases when a > 2 follow inductively by using the identity: 

A^[A- a - 2 ,ZW = A /3 - 2 [A- a ,Z]A 7 + A /3 [A- 2 ,Z]A 7 " Q 
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The cases when a < follow by similar identities The proof of Proposition B.2 is complete. 

Proof of Proposition 4.3. The Proposition 4.3 is a simple consequence of Proposition B.2. 
Since Z is in y , the claim 2) is covered by the bound on (B.3). We next prove 3). The operator 
K is in y 2 and Z is in y l . Power counting would suggest that [K , Z] e y 2 . However, by 
Condition HI, we know that V F equals a linear term plus a term in y°, and hence explicit 

computation shows that [K Q , Z] <E y l . Hence the assertion is covered by the bound on (B.3). 
Since M e y° we see from Lemma B.l that [M, A 2 ] e y 1 , and therefore the claim 1) follows 
from the bound on (B.4). Using again explicit calculation and Condition HI we see that [Z, A 2 ] 
is in y l (and not only in y 2 ) and [K Q , A 2 ] is in y 2 (and not only in 3^ 3 ), and hence the cases 4) 
and 5) follow by applying again the bound on (B.4). The proof of Proposition 4.3 is complete. 

Appendix C: A Generalized Core Theorem 

We prove here the following result from Section 4: 

Theorem 4.8. Let B be a Banach space. Let A : D(A) — > B be m-accretive. For all 
n = 1, 2, . . ., if V is a subset of D(A n ) and is dense in B and furthermore V is invariant under 
the semi-group e At , then V is a core for A n . 

Proof. Let |||/||| n = ^=o \\A* f\\. Then one has 

(i) (D(A n ), HI • IIIJ is complete. 

(ii) e At is a strongly continuous semi-group on (D(A n ), \\\ ■ ||| n ). 

The statement of the theorem is equivalent to the following: If V denotes the closure of 
V in the norm ||| ■ ||| n , then we have 

V n = D(A n ) . (C.l) 

We show this by induction. For n = 1, this is the core Theorem, [Da], Thm 1.9. Let us assume 
that (C.l) holds for n — 1. Let / e D(A n ), so there is a sequence {f m } £ V such that 

lim |||/ m -/||| B _, =0. 

m^oo 

With the f m we construct a sequence which converges to / in the ||| • ||| n norm. We set 

J m 

By the above property (ii), e As is strongly continuous in s in D{A n ) and hence G V . 

We set 

f(n,t) = [\ (t~s) n ~ l A 

1 Jo (n-1)! f - 



Jo (n-1)! tm ' 
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Since e As g is n-times differentiable in s when g e D(A n ), we obtain, upon integrating by parts: 



?(n,t) _ J?(n,t) 



n-1 

In ~~ ^ ] 

j=0 



+ 



o (n-1)! 



e^(/ m -/) 



+ 



(n-1) 

E 

i=o 

t n-l 



■M- '(/,„-/) II 



,n /("-I) 
< ^-111^-/111^ + ^ + 1)^-^11^- '(.,,„ 

= o(l) as m — > oo , 



by the inductive hypothesis. This shows that /( n '*) e V . 

To conclude we show that \\\^f M - f\\\ n -> as t -> 0. We have 



™ f(n,t) 



n-1 . „ 

= E - / 



*dg (t s Tj- ( e As -l)A j f 



+ 



n! 



(n-1) 



-(/ d g v ; ^— e ^M"/-l-A n /) 
t-Vo (n-1)! ■■' ' 



n!" 



Using that e A */ is n-times differentiable in t, that 



o (n-1)! 



is the remainder term in the Taylor expansion of e At f, and that for g e -D(-A) 

(e As -l^= f S due Au Ag, 
Jo 

if g E D{A) we obtain the bound 



1^1 /(»,*) 



n-1 

= E 



™ / d, 



(t-8) 



n — l ps 



(n-1)! 



due Au A J+l f 



+ 



n! 



f-f-Af -^/) 



n! 



< 0(t) +o(l) . 



Anharmonic Chains Coupled to Two Heat Baths 



43 



This shows that / e V . This proves that V = D(A n ) as required. 
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